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Abstract: 

Since 2002 the VHE (>100 GeV) 7-ray flux of the high-frequency peaked BLLac PKS 2155-304 has 
been monitored with the High Energy Stereoscopic System (HESS). An extreme 7-ray outburst was de- 
tected in the early hours of July 28, 2006 (MJD 53944). The average flux above 200 GeV observed during 
this outburst is —7 times the flux observed from the Crab Nebula above the same threshold. Peak fluxes 
are measured with one-minute time scale resolution at more than twice this average value. Variability is 
seen up to —600 s in the Fourier power spectrum, and well-resolved bursts varying on time scales of —200 
seconds are observed. There are no strong indications for spectral variability within the data. Assuming 
the emission region has a size comparable to the Schwarzschild radius of a — 10 9 M@ black hole, Doppler 
factors greater than 100 are required to accommodate the observed variability time scales. 



Introduction 

In the Southern Hemisphere, PKS 2155-304 (red- 
shift z — 0.116) is generally the brightest blazar at 
VHE energies, and is probably the best-studied at 
all wavelengths. The VHE flux observed [1] from 
PKS 2155-304 is typically of the order -15% of 
the Crab Nebula flux above 200 GeV. The highest 
flux previously measured in one night is approx- 
imately four times this value and clear VHE-flux 
variability has been observed on daily time scales. 
The most rapid flux variability measured for this 
source is 25min [2], occurring at X-ray energies. 
The fastest variation published from any blazar, at 
any wavelength, is an event lasting —800 s where 
the X-ray flux from Mkn501 varied by 30% [16] 1 , 
while at VHE energies doubling time scales as fast 
as —15 minutes have been observed from Mkn 421 
[9]. 



As part of the normal HESS observation program 
the flux from known VHE AGN is monitored reg- 
ularly to search for bright flares. During the July 
2006 dark period, the average VHE flux observed 
by HESS from PKS 2155-304 was more than ten 
times its typical value. In particular, an extremely 
bright flare of PKS 2155—304 was observed in the 
early hours of July 28, 2006 (MJD 53944). This 
contribution focuses solely on this particular flare, 
which is described in more detail in [4] . 

Results from MJD 53944 

A total of three observation runs (—28 min each) 
were taken on PKS 2155—304 in the early hours 
of MJD 53944. These data entirely pass the stan- 



1. Xue & Cui [16] also demonstrate that a 60% X-ray 
flux increase in —200 s observed [8] from Mkn 501 is 
likely an artifact. 
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Figure 1: I(>200 GeV) observed from 
PKS 2155— 304 binned in 1-minute intervals. 
The horizontal line represents I(>200 GeV) 
observed [3] from the Crab Nebula. The curve 
is the fit to these data of the superposition of five 
bursts (see text) and a constant flux. 
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Figure 2: The Fourier power spectrum of the light 
curve and associated measurement error. The grey 
shaded area corresponds to the 90% confidence in- 
terval for a light curve with a power-law Fourier 
spectrum P v oc v~ 2 . The horizontal line is the av- 
erage noise level. 



dard HESS data-quality selection criteria, yielding 
an exposure of 1 .32 h live time at a mean zenith an- 
gle of 13°. The analysis method is described in [4]. 
The observed excess is 1 1771 events (168cr), corre- 
sponding to a rate of ^2.5 Hz. This is the first time 
the detected VHE 7-ray rate has exceeded 1 Hz. 

Flux Variability 

The average integral flux above 200 GeV ob- 
served from PKS 2155-304 is I(>200 GeV) 
= (1.72±0.05 sta t±0.34 syst ) x 10- 9 cm- 2 s-\ 
equivalent to ~7 times the I(>200 GeV) observed 
from the Crab Nebula (Icrab, [3]). Figure 1 shows 
I(>200 GeV), binned in one-minute intervals, 
versus time. The fluxes in this light curve range 
from 0.65 Icrab to 15.1 Icrab, and their fractional 
root mean square (rms) variability amplitude [15] 
isF var = 0.58 ±0.03. This is ~2 times higher than 
archival X-ray variability [18, 19]. The Fourier 
power spectrum calculated from Figure 1 is shown 
in Figure 2. There is power significantly above 
the measurement noise level up to 1.6 x 10 -3 Hz 
(600 s). The power spectrum derived from the 
data is compatible with a light curve generated 
by a stochastic process with a power-law Fourier 
spectrum of index -2. An index of -1 produces too 
much power at high frequencies and is rejected. 



These power spectra are remarkably similar to 
those derived in X-rays [18] from the same source. 

Figure 1 clearly contains substructures with even 
shorter rise and decay time scales than found in 
the Fourier analysis. Therefore, the light curve 
is considered as consisting of a series of bursts, 
which is common for AGN and 7-ray bursts 
(GRBs). To characterize these bursts, the "gen- 
eralized Gaussian" shape from Norris et al. [13] 
is used, where the burst intensity is described by: 
I(i) = Acxp[-(|t - £ max |/cr r: d) K ], where i max 
is the time of the burst's maximum intensity (A); 
cr r and <7d are the rise it < i max ) and decay 
(t > i max ) time constants, respectively; and n is 
a measure of the burst's sharpness. The rise and 
decay times, from half to maximum amplitude, are 
T r ,d = [ln2] 1 / K CT r ^- Five significant bursts were 
found with a peak finding tool based on a Markov 
chain algorithm [12]. The data are well fit 2 by a 
function consisting of a superposition of an iden- 
tical number of bursts plus a constant signal. The 
best fit has a \ 2 probability of 20% and the fit pa- 
rameters are shown in Table 1 . Interestingly, there 
is a marginal trend for k to increase with subse- 
quent bursts, making them less sharp, as the flare 
progresses, which could imply the bursts are not 
stochastic. The k values are close to the bulk of 

2. All parameters are left free in the fit. 
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Table 1: The results of the best x 2 fit of the su- 
perposition of five bursts and a constant to the data 
shown in Figure 1 . 
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41.0 


173±28 


610±129 


1.07±0.20 


58.8 


116±53 


178±146 


1.43±0.83 


71.3 


404±219 


269±158 


1.59±0.42 


79.5 


178±55 


657±268 


2.01±0.87 


88.3 


67±44 


620±75 


2.44±0.41 



those found by Norris et al. [13], but the time scales 
measured here are two orders of magnitude larger. 

During both the first two bursts there is clear dou- 
bling of the flux within r r . Such doubling is some- 
times used as a characteristic time scale of flux 
variability. For compatibility with such estimators, 
the definition of doubling time, T 2 = |LyAT/AI|, 
from [18] is also used 3 . Here, AT = Tj - T t , 
AI = lj - T, Ly = (Ij + Ii)/2, with T and I being 
the time and flux, respectively, of any pair of points 
in the light curve. The fastest T 2 = 224 ± 60 s is 
compatible with the fastest significant time scale 
found by the Fourier transform. Averaging the five 
lowest T 2 values yields 330 ± 40 s. 

Spectral Analysis 

Figure 3 shows the time-averaged photon spectrum 
for these data. The data are well fit, x 2 = 17.1 
for 13 degrees of freedom (d.o.f.), by a broken 
power-law function with £ B =430±22±80 GeV, 
ri = 2.71±0.06±0.10, and T 2 = 3.53±0.05±0.10. 
For each parameter, the two uncertainties are 
the statistical and systematic values, respec- 
tively. The time-averaged spectrum (T = 3.32) of 
PKS 2155-304 measured in 2003 [1], multiplied 
by the ratio (48.7) of I(>200 GeV) from the re- 
spective data sets, is also shown in Figure 3. De- 
spite a factor of ~50 change in flux there is qual- 
itatively little difference between the two spec- 
tra which is surprising. The lack of any strong 
(AT > 0.2) temporal variability in the VHE spec- 
trum within these data (tested on time scales of 5, 
10 and 28 minutes) is also surprising. 
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Figure 3: The time-averaged spectrum observed 
from PKS 2155-304 on MJD 53944. The dashed 
line is the best x 2 fit of a broken power law to the 
data. The solid line represents the fit to the time- 
averaged spectrum of PKS 2155-304 from 2003 
[1] scaled by 48.7. 

Discussion 

It is very likely that the electromagnetic emission 
in blazars is generated in jets that are beamed and 
Doppler-boosted toward the observer. Superlu- 
minal expansions observed with VLBI [14] pro- 
vide evidence for moderate Doppler boosting in 
PKS 2155—304. Causality implies that 7-ray vari- 
ability on a time scale t var , with a Doppler factor 4 
(S), is related to the radius (R) of the emission zone 
by R < ct val S/(l + z). Conservatively using the 
best-determined rise time (i.e. r r with the smallest 
error) from Table 1 for i var = 173 ± 28 s limits the 
size of the emission region to RS -1 < 4.65 x 10 12 
cm < 0.31 AU. 

The jets of blazars are believed to be powered by 
accretion onto a supermassive black hole (SMBH). 
Thus accretion/ejection properties are usually pre- 
sumed to scale with the Schwarzschild radius i?s 
of the SMBH, where R s = 2GM/c 2 , which is 
the smallest, most-natural size of the system (see, 
e-g-, [7])- Expressing the size R of the 7-ray emit- 
ting region in terms of Rs, the variability time 

3. Only values of T2 with less than 30% uncertainty 
are considered. 

4. With 8 denned in the standard way as [F(l — 
13 cosf?)] -1 , where T is the bulk Lorentz factor of the 
plasma in the jet, f) = v/c, and 6 is the angle to the line 
of sight. 
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scale limits its mass by M < (c 3 t val 5/2G(l + 
z))R s /R ~ 1.6 x 10 7 MqSRs/R. The re- 
ported host galaxy luminosity Mr = —24.4 (Ta- 
ble 3 in [11]) would imply a SMBH mass of or- 
der l-2xl0 9 M Q [6], and therefore, <5 > 60 - 
120i?/i?s- Emission regions of only a few R$ 
would require values of 5 much greater than those 
typically derived for blazars (<5~10) and come 
close to those used for GRBs, which would be a 
challenge to understand. 

Although the choice of a ^3 minute variability 
time scale in this article is conservative, it is still 
the fastest ever seen in a blazar, at any wavelength, 
and is almost an order of magnitude smaller than 
previously observed from this object. The vari- 
ability is a factor of five times faster than VHE 
variability previously measured from Mkn 421 [9] 
and comparable to that reported from Mkn 501 [5]. 
However, in terms of the light-crossing time of 
the Schwarzschild radius, R$ / c, the variability of 
PKS 2155-304 is more constraining by another 
factor 5 of « 6 - 12 for Mkn 421, and a factor of 
w 2.5 -5 for Mkn 501. 

The light curve presented here is strongly over- 
sampled, allowing for the first time in the VHE 
regime a detailed statistical analysis of a flare, 
which shows remarkable similarity to other longer 
duration events at X-ray energies. More detailed 
discussion of this outburst can be found in [4], 
and the event continues to be investigated with 
other statistical techniques. As the sensitivity of 
VHE instruments continues to improve, it is likely 
that similar extreme flaring episodes will be more 
commonly detected in the future. Similar flares 
will only strengthen the conclusion that either very 
large Doppler factors can be present in AGN jets, 
or that the observed variability is not connected to 
the central black hole. 
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